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Abstract 

In resent years, Vicsek model has attracted more and more attention and been 
well developed. However, the in-depth analysis on the convergence time are scarce 
thus far. In this paper, we study some certain factors that mainly govern the con- 
vergence time of Vicsek model. By extensively numerical simulations, we find the 
convergence time scales in a power law with In N in the noise-free case, where r 
and N are horizon radius and the number of particles. Furthermore, to accelerate 
the convergence, we propose a new model in which the speed of each particle is 
variable. The convergence time can be remarkably shortened compared with the 
standard Vicsek model. 
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1 Introduction 



In nature, collective motion of large numbers of organisms is one of the most 
ubiquitous biological phenomena, from motion of groups of ant [1], colonies of 
bacteria [2] and cells [3] in a small scale, to migration of flocks of birds and 
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schools of fish [4] in a large scale. Those different forms of collective behaviors 
can only be understood by considering the very number of large interactions 
among group members. Studies on this issue is significant. On one hand, we can 
deepen our understanding of such collective behaviors; on the other hand, we 
could probably extract some generic rules from those natural systems, which 
can be applied to other relevant realms, such as the control of unmanned 
vehicles or robots [5]. 

Inspired by biological collective motion, Vicsek et al. [6] described each in- 
dividual in the collective motion as an self-driven particle moving with a con- 
stant speed and adjusting the direction according to the average direction 
over neighborhood. By this means, each particle will achieve the same veloc- 
ity through finite steps [7] , namely a kinetic phase transition from disorder to 
order state [6,8,9]. Besides, many modified models about self-driven swarm are 
also proposed. For example, Couzin et al. [10,11] proposed a more biologically 
realistic model, and studied the effect of the repulsion, alignment, attraction 
and leadership on swarm formation. 

In a word, though its simplicity, the Vicsek model is of great academic sig- 
nificance [12], which has become a theoretical templet for the consensus of 
swarm. However, previous works mostly focus on the analysis of steady state, 
while they pay little attention on the time used to achieve that steady state. 
As a new aspect of investigating the collective motion, the convergence time 
also has great significance to be explored. On one hand, it can be a criterion 
of evaluation of a consensus strategy. On the other hand, detailed discussions 
on the relationship between the factors and the convergence time can provide 
us some insights about the collective dynamics. Therefore, in this paper, we 
will study which factors influence the convergence time and how can a group 
of agents get to be coherent as quickly as possible. First, we will discuss the 
relationship between the convergence time and the particle density as well as 
the horizon radius of individuals. Then, considering the fact that the speed 
of an individual in biological group or removable robots usually can vary in a 
certain range, we abandon the strong assumption in the Vicsek model, namely, 
the constant speed. Instead, we propose a swarm model with variable speed, 
upon which a new communication protocol is designed, and the convergence 
time under this protocol is shorter than that of the standard Vicsek model. 



2 Convergence time in the Vicsek model 

A group of N particles are considered which are moving in an L x L square 
with periodic boundary conditions. In the Vicsek model, the particles are 
moving in an identical constant speed but different directions, with the initial 
conditions that the particles are randomly distributed in the square, and their 
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initial moving directions are uniformly distributed in the interval [— tt, tt). At 
each time step, the direction of each particle is determined by the average 
direction of the velocity of all the particles (including itself) within the circle 
of horizon radius r centered the given particle. Mathematically speaking, the 
position of the ith particle is updating according to: 

Xi{t + St) = Xi{t) + Vi{t)dt, (1) 

while the corresponding discrete presentation is: 

Xi{t + l)^Xi{t)+Vi{t). (2) 

And its direction is updating as: 

e,{t + i) = {e,{t))r + Aei, (3) 

where A9i denotes the noise, and {9i{t))r denotes the average direction of all 
the particles within the horizon radius r, including the ith particle itself. A^j is 
a random variable uniformly distributed in the interval [—r]/2, r)/2\. Obviously, 
{6i{t))r is given by: 

tanmm = ^p^^. (4) 

{ViC0s9i{t))r 



Moreover, in order to measure the degree of consensus of all the particles, 
an order parameter [6] is introduced as 
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0<(j)<l. 



(5) 



A larger value of (f) indicates a better consensus, especially when = 1, all par- 
ticles move towards the same direction. In the circumstances of high density 
and low noise, all the particles will definitely approach the consensus state, 
namely reach the same direction of velocity after finite time steps [7] . 

Because of the limitation of horizon radius, each particle can only communi- 
cate with partial particles within the radius and change its direction according 
to this local information. Different horizon radii and particle densities will give 
rise to diverse convergence behaviors and times, which will be systematically 
investigated later. Without lose of generality, the area is fixed, thus the parti- 
cle density can be directly represented by the number of particles, A^. 

The simulation results about the convergence time T based on the noise- 
free Vicsek model {f]=0) are shown in Fig. 1, where T represents the number 
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Fig. 1. (color online) Dependence of convergence time T on horizon radius r and the 
number of particles N. In the simulation, all the particles move in a square shaped 
plane of linear size L = 5, with a constant speed as v = 0.05 . Control parameters r 
and A^, respectively, vary from 0.5 to 1 and from 100 to 500. The convergence time 
is obtained from the average over 500 independent runs. 



of time steps taken before the order parameter (f) firstly reaches the thresh- 
old (pt = 0.95. As long as the threshold (pt is larger enough (close to 1), the 
variance of its specific value will not impact much on the qualitative results 
shown in this paper. Fig. 1 clearly indicates that the convergence time de- 
creases with the increase of radius and density. This fact can be explained in 
the following two aspects, (a) Given the density of particles, the larger the 
horizon radius, the less steps taken to reach synchronization, because at each 
time step each particle can receive more information from others and thus 
make the adjustment of velocity more comprehensively, namely more close to 
the final convergence velocity, than those with shorter radius, (b) Given hori- 
zon radius, when the number of particles increases, although the percentage of 
particles which communicate with a given particle does not increase as well, 
higher density is also helpful to reduce the convergence time, because the par- 
ticle is making more integral adjustments at each time step for the increase of 
the number of other particles inside its horizon (those particles possess differ- 
ent directions of velocity, come from different time step before, bring 
in information from different areas at latest time step and then pass the in- 
formation of current adjusted velocity to different areas at the following time 
step). 
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Fig. 2. Convergence time T as a function of InN. The data points can be well 
fitted linearly in double logarithmic coordinates, with slope and error bars marked 
in the plot. Each data point is collected from the average over 500 independent runs. 



1200 



►-1000 

0) 

E 

■■f= 800 



o 

I 600 

CD 

I 400 
o 

200 

0-1 




1 •- 



-I 1 1 1 1- 



0.2 0.4 0.6 0.8 1.0 1.2 



Fig. 3. Convergence time T as a function of the noise rj. We set r =1.0 and iV=100. 

Through extensive simulations we eventually come to a numerical conclusion 
that the convergence time T follows a power law with In N. Fig. 2 shows a 
numerically approximate relation: 



r~ (r^lnAT) 



-1.29 



(6) 
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In the case with noise, namely rj > 0, the convergence time will increase with 
the noise strength. In Fig. 3. when r] < 1.1, the convergence time increases 
slightly, while when rj > 1.1, the convergence time increases dramatically, 
especially when rj > 1.2, the order parameter can not achieve 0t=O.95, namely 
with high noise this collective motion can not achieve consensus. This result 
is quite in according with Fig. 2 in Ref. [6]. Though quantitatively different, 
the sharp increasing of T (as the increasing of 0) can be generally observed 
for different parameters {r,N,(j)t). 



3 Fast convergence collective motion with variable speed 



Without taking account of generation and annihilation of particles, the num- 
ber of particles would be generally fixed. In addition, it might be difficult to 
accelerate the consensus process by simply changing the horizon radius, for 
which means higher requirement on both technology and the hardware cost. 
Considering the speed of particle in collective motion systems can be variable 
[13,14], we propose a new consensus strategy that achieves improvement in 
terms of the convergence time compared with the standard Vicsek model. 

With limited horizon radius, a particle could only make the judgment by 
itself from the local information it receives. In a completely chaotic case, al- 
though each particle updates its direction to the average one in its local area, 
this average direction may be far different from the final synchronization direc- 
tion. Therefore, it should be cautious to act with a comparatively conservative 
strategy, that is, taking a relatively lower velocity to prevent from the unnec- 
essary change of its position. Here the unnecessary change means that, if a 
particle changes its position hastily in order to communicate with another 
group of particles under such chaotic circumstance, it would be sure that such 
impatient behaviors are of no help. On the contrary, if all the particles choose 
the conservative strategy of moving with a comparatively lower velocity, then 
every particle could have sufficient time to communicate with its neighbors, 
leading to a possibly faster convergence. Only when a certain moving direc- 
tion is dominant among its neighbors of a particle, it can surely update its 
direction as that one and take a relatively higher velocity since in that case it 
is unnecessary for this particle to continue hesitating. Therefore, a particle's 
velocity should be somehow determined by the degree of its local consensus. 

In order to present the degree of local synchronization quantitatively, 0j is 
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introduced, called the local order parameter of consensus for the ith particle: 
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, O<0, <1, 



(7) 



where Ni is the number of particles within the horizon radius of the ith particle 
(including itself). The larger value of the higher degree of local consensus 
among the neighbors of the ith particle. Especially when (f)i{t) — 1, all i's 
neighbors move towards the same direction. 

Besides, in order to compare with the Vicsek model in terms of the con- 
vergence time, we suppose that the absolute speed can vary from to 0.05. 
According to the discussion mentioned above, this consensus strategy with 
variable speed should satisfy: 

a) When all particles in i's neighborhood arrive at an ordered direction as 
(l)i{t) = 1, Vi{t + 1)^0.05; 

b) When — 0, namely the motions of the particles in i's neighborhood 
is completely disorder, Vi{t + 1) approaches to zero. 



Consequently, we set the speed of the ith particle as: 



(8) 



where Vmax is set as 0.05 in this paper. Here /3 is a free parameter. When P — 
the protocol degenerates to the standard Vicsek model, while for > 0, a 
particle will move faster in a better synchronized local circumstance. Indeed, 
in the present protocol, speed not only determines the position in the next 
time step, but also serves as a carrier, possessing the information about the 
local order parameter. The moving direction of the ith particle is also updated 
following Eq. (4). Note that, when approaches to 1, 0j gets close to 1 as well. 
Therefore from Eq. (8) we know that the speeds of all particles at this time 
are all close to 0.05, which automatically ensure consensus on the absolute 
speed. 

To sum up, the present protocol with adaptive speed can be described as 
follows: 

(1) Determine the initial position and speed of every particle i; 

(2) Evaluate the local order parameter of each particle, and determine its 
next direction and speed according to Eqs. (4) and (8); 

(3) Calculate the current order parameter of all particles; 

(4) Repeat (2) and (3) until (p approaches to 1. 
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Fig. 4. (color online) Convergence time T as a function of In N with different /5. 
The data points can be well fitted linearly in a double logarithmic coordinates, with 
slopes and error bars marked in the plot. Control parameters r and N, respectively, 
vary from 0.5 to 1 and from 100 to 300. The convergence time is obtained from the 
average over 500 independent runs. 

In the numerical simulations, we first consider the situation without noise 
[rj = 0), and still assume that all particles move in a square-shaped plane 
of linear size L = 5. The relation between the convergence time T and the 
horizon radius r and the number of particles N with different (3 is shown in 
Fig. 4. From this figure, we can see that in this new protocol the convergence 
time T also has a power function with In N, the same rule as illustrated in 
section 2. Especially, when (3 is not large enough {(3 < 5), the exponents are 
almost the same, nearly equal to 1.30. More importantly, in Fig. 4 we can find 
that, compared with the standard Vicsek model, the convergence time can be 
shortened, which demonstrate the advantage of this modified model in terms 
of the convergence time. 

To further explore this point, we keep the horizon radius r and the number 
of particles fixed and try to find out the optimal j3 subject to the shortest 
convergence time. In Fig. 5, we can see that in all the four curves the optimal 
/3 is 2 and this optimal /3 is not susceptive to the horizon r and the number 
of particles A^. Especially, the convergence time can be highly shortened when 
the horizon radius r is small. 

In Fig. 6, we can see that when (3 is small (e.g. f3 = 0.1), no matter what is the 
current speed distribution, the speeds of most particles in the next time will 
approach to the upper limit 0.05, resulting in almost the same performance 
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Fig. 5. (color online) Convergence time VS. (3 with different r and A^. 
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as the standard Vicsek model; besides, when /? is too large (e.g. /? = 10), the 
speeds of particles will tend to be 0, thus even when the neighbors of a given 
particle is relatively ordered, it still keeps rest and therefore cannot intercom- 
municate with others, thus resulting in a extremely long convergence time. 
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For a proper value of (5 (e.g. (3 = 2), when a particle has a large local order 
parameter, it will have a great confidence on its current direction, thus leaving 
the area with a high speed and using the information (speed direction) to in- 
fluence others; Otherwise, it will cast doubt on whether the current direction 
is up to the demand of consensus, and in order to avoid misguiding others 
with its unconfirmed direction, it will march on with a relatively slow speed. 



4 Conclusions 

The collective dynamics of intelligent agents is not only the extensive phe- 
nomena in nature, but also an important problem required in-depth investiga- 
tion in engineering. Most of the previous studies concentrated on the depiction 
and modeling of the swarm itself. The systematical analysis about the conver- 
gence time were rarely reported. In this paper we have studied the relationship 
between the convergence time and the particle density as well as the horizon 
radius in the Vicsek model, and found that the convergence time T has a power 
function of In N. Furthermore, we have designed a motion protocol variable 
speed, in which the speed is not only a moving parameter, but also can deliver 
some information about its situation of local synchronization. Under such pro- 
tocol, the convergence process can be remarkably accelerated compared with 
the standard Vicsek model. 
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